Magnetic dephasing in a mesoscopic spin glass 
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We study the conductance of a mesoscopic metallic conductor with magnetic impurities frozen in a 
spin glass phase. We describe the correlation of conductance between two frozen spin configurations 
in terms of dephasing rates for the usual low energy modes of weak localization theory. They 
depend on the overlap between the spin configurations and are obtained through the crossover 
between orthogonal and unitary universality classes of a nonlinear sigma model. 
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In a metallic conductor of /im size, the interplay be- 
tween quantum coherence of electrons and disorder is 
at the origin of several remarkable phenomena at low 
temperature [H, [3]. The interferences between differ- 
ent diffusive paths in the sample lead to a strong de- 
pendance of the conductance on the disorder realization 
through the geometry of these diffusive paths. In partic- 
ular the conductance exhibits universal sample to sam- 
ple fluctuations and reproducible random variations as a 
function of a transverse magnetic flux. These coherent 
transport phenomena are naturally suppressed by var- 
ious perturbations which are usually referred to as de- 
phasing sources. This includes inelastic scattering, for 
example from phonons, free magnetic impurities, other 
electrons. In this case two electrons following the same 
path encounter different potentials. This leads to a ran- 
dom relative phase between them and the suppression of 
interference effects. Such a dephasing source acts on the 
electrons themselves, and is usually taken into account 
phenomenologically through a dephasing rate 7^ for the 
electrons. Elastic scattering by symmetry breaking po- 
tentials by e.g random spin-orbit interactions or frozen 
magnetic impurities is also considered as another source 
of dephasing. In this case, the origin of dephasing is phys- 
ically different from that in the inelastic case. Quantum 
corrections to transport result from interference between 
electrons travelling along loops of diffusive paths either 
in the same directions (Diffuson modes) or in the op- 
posite directions (Cooperon modes). In the presence of 
magnetic disorder, the spins of counter-propagating elec- 
trons experience different sequences of rotations. As a 
result the corresponding interferences are gradually sup- 
pressed for longer and longer loops. This suppression 
is usually interpreted as dephasing of the correspond- 
ing Cooperon modes. Moreover the magnetic disorder 
selects out the Diffuson modes not affected by this de- 
phasing: the relative phase of two electrons with spins 
forming a singlet state is insensitive to this spin rotation 
sequence, in contrast to that of triplet states. Hence, 
this symmetry-breaking elastic scattering induces a spa- 
tial decay of these particular diffusion modes, effectively 
accounting for this electron's dephasing phenomenon. 

This magnetic dephasing of electrons and its signa- 



ture on coherent transport has been recently proposed 
as a promising probe of spin glass physics [3|. In the 
spin glass phase, which is a fascinating but poorly under- 
stood state of matter impurity spins {Si} coupled by 
frustrating interactions freeze below the transition tem- 
perature Tg. Of particular interest for its understanding 
are the correlations between different configurations of 
spins {^l^^} and {-S'f^^} corresponding to either different 
times (to probe aging) 0, Q, or different quenches be- 
low Tg to probe possible mean-field complexity. In this 
last case, the complexity of phase space is encoded in the 
different overlaps 912 (or distances) between spin config- 
urations In a given sample, a magnetoconductance 
trace carries a fingerprint of the disorder configuration, 
including the configuration of frozen spins {S^"'}- Fol- 
lowing the above discussion its properties depend on the 
magiretic dephasing in this spin configuration. Moreover, 
the correlation between conductances for different spin 
states (different quenches or times) is determined by the 
correlation of magnetic dephasing in the two correspond- 
ing spin states. Indeed for different spin configurations, 
two electrons prepared in a Singlet state, and travelling 
along the same path now experience two different mag- 
netic potentials. Hence the Diffuson Singlet mode built 
out of these different impurity spin configurations is now 
dephased. Moreover this dephasing is a measure of the 
distance or overlap between the spin configurations: the 
larger the number of spins differing in the two configu- 
rations, the stronger the relative dephasing of electrons. 
Hence, such measurement of correlations of conductance 
open the route for measurement of correlations between 
frozen spin configurations [3, 0| . This naturally requires 
a quantitative description of this relative magnetic de- 
phasing between different spin configurations, on which 
we focus in this Letter. 

Beyond the perturbative diagrammatic theory [2], a 
natural framework for describing the quantum correc- 
tions to conductance is the field theoretical nonlinear 
sigma model We will show that this method allows 
for a very efficient and elegant description of the magnetic 
dephasing, by treating all elastic scattering potentials on 
the same footing. Altland showed how the inelastic de- 
phasing effects can be phenomenologically accounted for 
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within this nonhnear sigma model ll|. While his ap- 
proach required the introduction of a new fictitious scalar 
potential, in the present context the source of dephasing 
is already an elastic scattering potential. Moreover, the 
dephasing accounts for the crossover behavior between 
different universality classes characterizing both the An- 
derson transition in systems of various statistical sym- 
metries and their universal weak localization properties 
fl^ . In particular they encode the number of indepen- 
dent diffusive modes contributing to the quantum cor- 
rection of conductance for weak disorder. In the present 
case, without magnetic impurities the statistics of con- 
ductance is described by the orthogonal class with degen- 
erate diffusive states. When adding magnetic disorder, 
the level degeneracy is lifted and simultaneously the uni- 
versality class is restricted to the unitary class. Studying 
the magnetic dephasing induced by frozen magnetic im- 
purities amounts to study the cross-over between these 
two universality classes. This is achieved by considering 
all massless modes of the orthogonal case, and describ- 
ing the different gap openings with increasing magnetic 
disorder, thereby extending non-perturbatively the work 
of pJi] to the spin glass physics of interest here. 

We consider a d dimensional mesoscopic metallic sam- 
ple of size L containing impurities inducing two different 
types of random scattering of the conduction electrons: 
(i) nonmagnetic random scalar potential V{r) coupled to 
the fermionic electronic field ip{r) as V{r)ip{r)ip{r). This 
potential is assumed to be Gaussian with {V{r)) = and 
variance 



{V{r)V{r')) 



1 



S{r-r'), 



(1) 



where vq is the one-electron density of states and Ty - 
the corresponding elastic mean free time, (ii) a mag- 
netic disorder U (r) originating from a collection of frozen 
magnetic impurities Si with coupling to the electrons 
U{r) ■ ip(r)aip{r) with a being the Pauli's matrices. U{r) 
is a three dimensional field taken Gaussian with zero 
mean in the spin glass state. Here we focus on both the 
conductance fluctuations in a single spin configuration, 
and conductance correlations between different spin con- 
figurations corresponding to different spin glass states. 
For simplicity we consider correlations between different 
realizations of magnetic potential indexed by u which de- 
pends only on the overlap q^u' as 

U^Ar)Uu'-Ar') = ^^^^^^/(^ - ^'), (2) 



where Quu = 1 and r, is the elastic mean free time for 
scattering by magnetic impurities. 

Following Ref. [3] we define the covari- 
ant and contravariant bispinors fj = (Crj)'^ = 
(— '01, — V^ij V'ii "V't) ' which are related by the 
charge-conjugation matrix C = iai (g) a2- For our 



purpose we need to introduce AN copies of the original 
system enumerated by: (i) p = R, A for generating the 
advanced and retarded Greens functions; (ii) different 
replica a — l,..,N for disorder averaging using the 
replica trick; (iii) u — 1,2 to compute the correlation 
of conductance between 2 different configurations of 
magnetic disorder. In terms of the bispinors the action 
for the noninteracting electrons can be written as 
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cfir fi{r) S,-V{r) - U{r) - T. r]{r), 



(3) 



where ^ = ep — p^/(2m) and — a. To com- 

pute the dimensionless conductance correlation in a given 
sample with two different spin configurations (AG')^^^, = 

{G{V, Uii)G{V, Uu'))c we generalize the formahsm of gen- 
erating functional introduced in Ref. 'IJI]. We define the 
current operator J — i / 2{ri ® dfj — drj ® ff) . Then the 
generating functional for the conductance reads Z[A\ = 
^ 'Dfj'Drj eyiY> \iS -\- /^Tr(JrAr)] (we omitted a demina- 
tor with is suppressed in the replica limit iV ^ 0). The 
symbol Tr stands for traces over all indices while tr will 
be used for traces only over retarded-advanced, replica 
indices and integration over r. The correlation of dimen- 
sionless conductance can be expressed as 
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(4) 

Here 51 = ctq T2 where ctq and T2 are respectively the 
2x2 unit matrix and matrix with all entries equal to 1. 

Averaging over disorder we obtain the quartic 
term S.m = l/(4i7ri^oTe) J^fjrjPfjr] with P^^{u,u') = 
{'''e/Tv)SjS'i + (re/3Ts)(7„„'E^E;, where we introduced the 
total mean free time which is not fixed at this stage. 
Introducing a new field Q of the same rank and symme- 
try as fj 1^1] we perform the Hubbard-Stratanovich trans- 
formation on the quartic term. The charge-conjugation 
symmetry automatically ensures decoupling in two pair- 
ing channels: (i) Diffuson modes corresponding in the 
standard diagrammatics to the ladder diagrams with a 
small transfer momentum, (ii) Cooperon modes corre- 
sponding to the ladder diagrams with a small sum mo- 
mentum 0. Integrating out the bispinor fields we ex- 
clude all " fast" modes and obtain the free energy in terms 
of "slow" degrees of freedom Q 



^QA-Q-^Trln 



,(5) 



where K satisfies P^iKli = '^"'^m- 

The first step in elucidating the physics of low energy 
excitations is to find out the classical solution correspond- 
ing to a spatially-homogenous field Q. The correspond- 
ing saddle-point condition ttvqKQ = + iQ/(2re)]~^ 
can easily be solved in the limit of no magnetic disorder 
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1/ts = 0. In this limit we replace QKQ TrQ^ so that 
the homogeneous solutions satisfies = 1 and TrQ = 0. 
The corresponding saddle-point manifold can be param- 
eterized by Q = (1 - A(l - where A = 
diag(l8Ar, — l8Af)p,p' and W is an anti-hermitian matrix 
anticommuting with A. For finite Ts, enforcing = 1 
and TrQ = for the classical solution leads to an addi- 
tional condition KQ = Q. This requirement controls the 
lowering of the symmetry from the orthogonal class to 
unitary one for 5„ti' — 1- Out of the 16 eigenvalues of the 
matrix K only two of them can be fixed to 1 by a proper 
choice of the elastic mean free time Tg. Through this 
procedure we recover the well known Matthiessen rule 
1/re = 1/Ty + 1/ts. Moreover, the corresponding eigen- 
modes Q of matrix K are diagonal matrices describing 
the Diffuson Singlet state, the only massless mode in the 
unitary ensemble. 

Having determined the classical solution for finite r^, 
we can now expand ([5]) around it. Performing a gradient 
expansion around the above homogeneous saddle point 
we rewrite the effective action as 



mO] = -^/Tr 



—WMW + D{dWy' 



(6) 



we introduced the matrix M'^-j.{u,u') — K^''i.\ 



k:{u,u) 



where D — VpTe/d is the usual diffusion constant and 

^, , K^l{u,u)] . The usual gauge symme- 
try [10| ensures that the dependence of ([5]) on enters 
only in the combination dQ — [A, Q] and can be easily 
restored in Eq. ([6]) . The matrices W are of the form 
W — offdiag(-B, — i3^)p^p' with B satisfying the charge- 
conjugation symmetry B = CB*C'^. To parametrize 
B, we introduce a generalization of the standard quater- 
nion basis to the bispinor space by defining the 
following basis fpif^+i, = c^^i^afj, (gi Ou where c = i for 
(0 < M < 2)0(1 < < 3) and for /i = 3, = 0, 
and 1 otherwise. These matrices i/)^ form a complete 
set and satisfy the charge conjugation symmetry and re- 



lation Tr[ 



~ 4(5^j,. In the resulting decomposition 



B = (pfj^b^ and i?^ — 



the 6p are now real matrices. 



Using this decomposition of matrices B, we obtain the 
following quadratic part of the free energy in terms of real 

variables 6^ : T[W, 0] = itvo E {Tm + Dq^} K^pb'^^f}- In 
deriving this expression, we used that the field 0^ natu- 
rally diagonalize the matrix A4, and satisfy the identity 
(j)jj^A4(j)u = 4Te(5p,y7^. This diagonalization provides the 
masses (or dephasing rates) 7^ of the diffusing modes 
beyond the previous first order expansion in Ty/rg Q. 
As expected, among all eigenvalues 7^ there are only 4 
different masses 



(AG) 




FIG. 1: The conductance fluctuations (G^) — (G) for diffusive 
wire (ID), as a function of the wire's length L/Ls, and for 
different values of relative amplitude t^/ts of magnetic versus 
scalar disorder. The dashed curve corresponds to perturbative 
dephasing results of These function extrapolate between 
the usual universal values 8/15 (orthogonal class) and 1/15 
(unitary class). The inset shows the contributions of the 4 
different Cooperon/Diffuson modes for r„/rs = 0.2. 
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FIG. 2: The conductance correlations {G{V, Si)G{V, S2)) - 
2 

{G{y,Si)) as a function of the overlap 1 — g between the 
different magnetic disorder (spin) configurations Ui and U2 
for a ID conductor. The relative strength of both disorders is 
t/ts = 0.3, and wire length L/ Ls — 5 (left) and L/Ls = 1.8 
(right). The relative contributions of Cooperon and Diffu- 
son Singlet and Triplets modes show that while all contribute 
significantly to the overall value of the correlations, their gi2 
dependance is dominated by the Diffuson Singlet term. 



7b = 7r, 



1 + KBTT^t- 



(7) 



where 7™ = Ty/{TeTs), = 1 



-1/3, K 



1/3. These masses can now be identified with 



7?, 



7?, 



7sS 



7^ for the Singlets 



the dephasing rates 
and Triplets components of usual Diffuson and Cooperon 
modes. Indeed, in a single disorder reahzation = 1), 
one recovers that longitudinal fluctuations in the classical 
manifold are massless : they correspond to Diffuson Sin- 
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glet mode. The expansion now also provides the masses 
for the modes transverse to this manifold : the Diffuson 
Triplets and Cooperon modes. These masses, or elastic 
dephasing rates, are valid as long as the above classical 
solution is stable, implying at least r^/r^ < 1 from the 
above expression. The expressions ([7]) also provide the 
dephasing between different magnetic disorder configura- 
tion {qu,u' 7^ 1), which are now all finite. 

Let us now turn to the study of the correlation of 
conductance. The necessary terms of order in the 
generating functional read to one loop order: Z[A] = 
^iTTiyoD)^[{FnFn)o + {Fi2Fi2)o + 2{FuFi2)Q] with ver- 
tices Fii = J^Tt[WAWA] and F12 = J^Tt[WWAA]. 
Expressing W in the basis of performing Wick's con- 
tractions of 6^ and substituting in Eq. ^ we obtain the 
correlation of conductance. The result can be conve- 
niently written in terms of the four dephasing lengths 

given by L^{u, u') = J D/{2j^) and = ^/DT^ : 



(AG)L' = (G(T/,{/„)G(T/, {/„')> - (G(T/,C/„)) 



f 



3/ 



+ / 



3/ 



(8) 



with fix) = iJ2q^o[iLqf +x^/2]-^. Specifying this 
result to the case d = 1 of a diffusive wire [ij], we use 
fix) = 3a;-''[x2csch2 (a;/V^)+V2a;coth (x/V2)-4]. For 
Quu' — 1 the scaling function ^ shown for d = 1 in Fig. [1] 
gives the sample to sample conductance fluctuations and 
describes the crossover between the orthogonal and uni- 
tary universality classes. The inset shows that all four 
terms contribute to the scaling function but the correc- 
tion to the unitary fixed point is dominated by Diffuson 
and Cooperon Triplets. In Fig. [2] we plot the average 
correlation between conductances for different magnetic 
disorder realizations, experimentally measurable through 
correlation of magnetoconductances [3]. In the typical 
regime of wire length L large compared to magnetic de- 
phasing lengths of order Ls, the correlations of conduc- 
tance decay as we lower the overlap between the corre- 
sponding spin configurations. Moreover, this decay as a 
function of I — g is dominated by the dephasing of the 
Diffuson Singlet contribution while other modes almost 
compensate for each other (see left side of Fig. [5]). Note 
that in the opposite regime L < L^, anomalous behavior 
can appear with a qi2 dependance of (AG)^2 dominated 
by the Diffuson Triplet and Cooperon Singlet contribu- 
tions leading to a (small) increase as a function of l — qi2- 
Finally, let us stress that this monotonous decrease 
of (AG)^^, in the experimental regime of interest 15 1 
L > Ls allows for an interesting and unique test of a 
spin glass mean-field theory. Indeed, this theory predicts 
the ultrametricity of the spin glass phase space in the 
thermodynamic limit. According to this prediction, if 



we consider three spin configurations (or 3 and sort 
their mutual overlap according to 1 — qi2 > 1 — 513 ^ 1 ~ 
(723, then 1 — ^12 — l^qis- This condition easily translate 
into the practical test (AG)?2 = (AG)?3 if (AG)?2 < 
{AG)l, < (AG)i3. 

To conclude, we have shown how to account natu- 
rally for the magnetic dephasing of diffusing electrons 
within the usual nonlinear sigma model. Motivated by 
the study of mesoscopic spin glass wires, we have used 
this formalism to study the relative dephasing rates be- 
tween different magnetic disorder configurations, and the 
corresponding correlations of conductance fluctuations 
amenable to direct experimental measures. Let us finally 
stress that an advantage of this field theoretical method 
is its flexibility, allowing for interesting extensions includ- 
ing the incorporation in our approach of more complex 
statistical correlations of spin configurations, along the 
lines of , as well as higher moments of the conductance 
correlations. 
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